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ABSTRACT 

We have identified a complete, flux-limited, (Sieo > 120 mJy), sample of 160/xm- 
selected sources from Spitzer observations of the 1-deg 2 ISO Deep Field region in the 
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Lockman Hole. Ground-based UV, optical and near-infrared (NIR) photometry and 
optical spectroscopy have been used to determine colors, redshifts and masses for the 
complete sample of 40 galaxies. Spitzer-IRAC+MIPS photometry, supplemented by 
ISOPHOT data at 90fim and 170/um, has been used to calculate accurate total infrared 
luminosities, Lir(8 — 1000/mi), and to determine the IR luminosity function (LF) of 
luminous infrared galaxies (LIRGs). The maximum observed redshift is z ~ 0.80 and the 
maximum total infrared luminosity is log (Lih/Lq) = 12.74. Over the luminosity range 
log (Liji/Lq) = 10 — 12, the LF for LIRGs in the Lockman Hole Deep Field is similar to 
that found previously for local sources at similar infrared luminosities. The mean host 
galaxy mass, log (M/M ) = 10.7, and dominance of HH-region spectral types, is also 
similar to what has been found for local LIRGs, suggesting that intense starbursts likely 
power the bulk of the infrared luminosity for sources in this range of LrR. However for 
the most luminous sources, log (Ljr/Lq) > 12.0, we find evidence for strong evolution 
in the LF cx (1 + z) 6±1 , assuming pure number density evolution. These ultraluminous 
infrared galaxies (ULIRGs) have a larger mean host mass, log (M/Mq) = 11.0, and 
exhibit disturbed morphologies consistent with strong-interactions/mergers, and they 
are also more likely to be characterized by starburst-AGN composite or AGN spectral 
types. 

Subject headings: galaxies: evolution, infrared: galaxies 

1. Introduction 

Deep surveys at rest-frame far-infrared (FIR) wavelengths are important for identifying sta- 
tistically complete samples of luminous infrared galaxies (LIRGs: Lir/L© > 11.0) - objects that 
appear to produce the bulk of the bolometric infrared luminosity density in the universe, and 
which are often "hidden" and misidentified in deep UV-optical surveys. Until recently, progress 
has been relatively slow in identifying complete samples of FIR sources selected at wavelengths 
A = 100 — 200/iin, which corresponds to the wavelength range where the majority of LIRGs at 
z < 1 emit their peak emission. The Infrared Astronomical Satellite (IRAS) all-sky survey pro- 
vided the first complete census of 60/mi-selected galaxies in the local universe (e.g. Soifer et al. 
1989; Sanders et al. 2003), but lacked the sensitivity and long wavelength coverage to detect LIRGs 
at z > 0.05. The Infrared Space Observatory (ISO) provided increased sensitivity and longer 
wavelength coverage, but background instabilities often limited the determination of exact source 
positions. The Spitzer Space Telescope eventually provided the combination of long wavelength 
sensitivity and background stability needed to detect sources at the ~ lOOmJy level with rela- 
tively accurate positions, and extensive follow-up observations are now underway to identify source 
properties and redshifts. 

In this paper, we report observations first begun as part of the U.S- Japan ISO-ISOPHOT 
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Deep Survey of a ~l-deg 2 region in the Lockman Hole (Kawara et al. 1998; Oyabu et al. 2005), 
and later expanded to include more recent infrared observations of the same region obtained as 
part of the Spitzer Wide-area InfraRed Extragalatic (SWIRE) survey (Lonsdale et al. 2003). Our 
final sample consists of MIPS-160/mi sources with fluxes greater than 120mJy. Multi-wavelength 
follow-up observations include Keck spectra of the majority of the sources, along with UV-optical 
photometry from the Sloan Digital Sky Survey (SDSS: York et al. 2000) for all of our sources, 
NIR photometry from the 2-Micron All Sky Survey (2MASS: Skrutskie et al. 2006) for most of our 
sources, and VLA 1.4GHz continuum images of the ISO-ISOPHOT deep fields (Yun et al. 2003). 
While other surveys (e.g. COSMOS: Scoville et al. 2007) may offer superior (particularly ancillary) 
data, the Lockman Hole survey discussed here subtends a total of 1.2 deg 2 and is thus comparable 
in total area to similar existing datasets. This survey therefore substantially increases the total 
area to date at these wavelengths and helps guard against cosmic variance. 

Sections 2 and 3 describe our field selection and identification of SWIRE MIPS- 160 sources, 
respectively. Methods used for identifying optical counterparts are presented in § 3, along with 
UV-NIR photometry and optical spectroscopy for individual sources. The spectral energy distri- 
butions (SEDs) and calculated infrared luminosities (Lir) for each source are presented in § 5, and 
these data are then used to determine the Infrared galaxy Luminosity Function (LF). Section 6 
presents evidence for evolution in the LF at the highest infrared luminosities, as well as a discussion 
of the properties of the host galaxies, including morphology, colors, masses and spectral types. Our 
conclusions are presented in § 7. 



2. Field Selection and FIR Observations 

The Lockman Hole (LH) is a region of the sky with the smallest HI column density (Lockman 
et al. 1986) and thus has been a favorite target for deep extragalactic surveys, particularly in the 
FIR where confusion noise from infrared cirrus is expected to be at a minimum. 

The LH was originally targeted for deep FIR observations with the Infrared Space Observatory 
(ISO) as part of the Japan/University of Hawaii (UH) ISO cosmology program that carried out 
observations of 2 small regions (LHEX and LHNW) using the ISOPHOT detector at 90^m and 
170/im. Both LHEX and LHNW are ~ 44' x 44' in size (see Figure 1) - the former is centered at 
a = 10 h 52 m 00 s , 5 = 57°00'00" and corresponds approximately to the field center of the ROSAT 
X-Ray Deep Survey Field (Hasinger et al. 1998), while the latter is centered at a = 10 h 33 m 55 s , 
5 = 57°46'20" and was chosen by the Japan/UH team to be centered on the true minimum HI 
column density position within the larger LH field. Note that a small sky area in the LHNW field 
was also mapped at 7^m using ISOCAM on ISO (Taniguchi et al. 1997) 

There is a fairly extensive published analysis of the ISOPHOT 90/um and 170/xm data for 
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both LHEX and LHNW, including initial source counts (Kawara et al. 1998), and final source 
lists (Kawara et al. 2004) , along with attempts to identify radio and near- infrared counterparts for 
individual sources (Yun et al. 2001; Sanders et al. 2001; Yun et al. 2003), and to obtain follow-up 
spectroscopy (Sanders et al. 2001; Oyabu et al. 2005). However, several factors prohibited using the 
ISOPHOT data for producing a well identified complete sample of FIR sources, first of which being 
the relatively large fluctuations in the ISOPHOT background, and the corresponding uncertainty 
in extracted source positions at both 90/im and 170/_,m. Further complications came with the 
realization that several radio and NIR sources were often found within the relatively large FIR 
beam. 

A major improvement in our ability to determine more reliable FIR fluxes and more accurate 
source positions was made possible once data from Spitzer were made public. In particular, SWIRE, 
a Spitzer Key Project (Lonsdale et al. 2003) included the LH as one of its deep survey areas 1 . 
SWIRE obtained maps in all 7 Spitzer IRAC+MIPS bands (3.6, 4.6, 5.8, 8.0, 24, 70, 160/mi), and 
covered both LHEX and LHNW, except for a small portion of the LHNW region missed by IRAC 
as shown in Figure 1. 

In this paper we use the SWIRE IRAC+MIPS data to produce a new flux-limited sample of 
FIR sources in the LHEX + LHNW fields. We also make use of all of our previous multi- wavelength 
imaging and spectroscopy along with new data from SDSS and 2MASS in order to first identify re- 
liable optical counterparts, and then to determine redshifts and total infrared luminosities for each 
source. We also use these data to characterize host galaxy properties (e.g. masses, morphologies, 
spectral types) in order to better understand the nature of the infrared galaxy population. 



3. The Lockman Hole 160^im Sample 

The MIPS-160^m detector samples emission at wavelengths closest to the rest- frame peak of 
the FIR/submillimeter SED of infrared luminous galaxies. Most galaxies detected at 160/iin are 
also detected at MIPS-70/um and MIPS-24/im as well as all four IRAC bands. For these reasons, as 
was the case with previous ISO observations of the LHEX and LHNW fields, we continue to focus 
on selecting a complete sample of extragalactic sources at the longest observed FIR wavelength. 

Figure 2 (left panel) shows the distribution of flux density of all ~500 sources in the SWIRE 
LH 160^m catalog. It is from this distribution that we take 120mJy as the flux limit for the sample. 
There are 39 SWIRE 160/mi sources within the LHEX and LHNW fields down to a flux limit of 
120mJy. The integral number counts from the catalog (Figure 2, right panel) follow a constant 



1 The SWIRE data release paper is available at: 
http:/ / swire. ipacxaltech.edu/swire/astronomers/publications/SWIRE2_doc_083105.pdf 
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slope at fluxes above the 120 mJy selection limit, and we therefore take the sample to be complete 
above this limit (e.g. Beichman et al. 1988; Soifer et al. 1989). 

We also compared the MIPS-160/iin catalog with our previous ISOPHOT 170/um catalog, 
and secondarily with the ISOPHOT 90/um catalog, and confirmed our previous suspicions that a 
significant fraction of the ISOPHOT 170/um sources were either spurious 3a noise peaks or weaker 
sources with boosted flux due to non-Gaussian noise fluctuations. However, we also discovered a 
few sources where our previous ISOPHOT data along with MIPS-70/^m data suggested that there 
should have been a MIPS-160^m source above 120mJy, yet there was no source listed in the MIPS- 
160/xm catalog. For these cases we went directly to the MIPS images and found that in each case 
there was a source in the 160//m image, so we extracted our own flux at 160/um by performing point- 
spread function (PSF) fitting on the image. This was done in consultation with the SWIRE team 
so that our flux extraction method would be in agreement with that used to construct the SWIRE 
catalog. In all but one case our extracted MIPS-160/um flux fell below the 120mJy completeness 
limit. However, for J105252. 76+570753. 7, the extracted MIPS-160^m flux of 149mJy was above 
the completeness limit, thus this source was added to the final list, bringing the total number of 
sources to 40. The distribution of these 40 sources across the LHNW and LHEX fields is shown in 
Figure 3. 

In addition, there are six sources listed in Table 1 that were included in the SWIRE MIPS- 
160/xm catalog but not the MIPS-70/um catalog. This was somewhat surprising given the sensitivity 
of the MIPS-70/im catalog, and the expected 70^m flux assuming even a fairly extreme 160/70 flux 
ratio. All six sources were in fact visible in the 70/um image, and hence, we again performed PSF 
fitting in consultation with the SWIRE team to extract a 70fim flux. Thus all 40 of our MIPS- 
160/im sources now have measured MIPS-70/um fluxes. 



4. Multiwavelength Data 

4.1. UV/Optical and NIR Images 

UV/optical and NIR imaging data are critical for identifying the sources responsible for the 
observed FIR emission, measuring redshifts and computing luminosities. The SWIRE data release 
actually includes optical (g'r'i') images obtained at Kitt Peak National Observatory (KPNO) cov- 
ering a large fraction of the LH field, with extracted photometry for those sources matched to the 
identified IRAC counterparts of MIPS-24//m sources (see below). However, SWIRE optical cover- 
age is uneven across our fields, a central portion of the LHEX (ROSAT) field was not targeted for 
followup and the southwest portion of the LHNW field falls beyond the edge of the IRAC survey 
(see Figure 1) and consequently was also not observed in optical. As a result, almost half (19/40) 
of our MIPS-160^m sources lack KPNO g'r'i' photometry in the SWIRE catalog. Fortunately, the 
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SDSS provides full coverage of the LH field, and we therefore make use of the ugriz catalog pho- 
tometry from Data Release 7 as well as display color-composite images for all of the MIPS-160/um 
sources. Additionally we use 2MASS catalogs to obtain JHK S photometry for the majority (27/40) 
of our MIPS-160/um sources. 

4.2. Identification of Counterparts 

As an initial step the 160/um sources were matched with their counterparts in the SWIRE 
70/xm, and IRAC+24 / um+Optical catalogs. A series of image cutouts from each of the Spitzer 
IRAC+MIPS bands plus optical images from KPNO were then assembled for each source in order 
to identify the correct optical counterpart (see Figure 4 for an example 2 ). The red, blue, and green 
circles overlaid on these images indicate the position of the centroid of each detection in the 160/im, 
70/um, and bandmerged IRAC+24/um catalogs, respectively. Each 160/tmi source has one and only 
one counterpart in the 70/um images, and nearly all (36/40) have only one counterpart at 24/im. 
In the cases with multiple 24/xm sources, we choose the source nearest the 160/um centroid. Table 
1 lists the source identification numbers from the SWIRE catalogs, as well as coordinates of each 
source taken from the SWIRE IRAC+24/xm+Optical catalog, except for sources lacking coverage 
in two or more IRAC bands (noted in the table) where we list coordinates from the SWIRE 24^m 
catalog. Following the identification of the Optical/IRAC counterpart for each of the MIPS-160/um 
sources, we then match these with objects in the 2MASS and SDSS catalogs and list their source 
IDs. In the case of 2MASS we prefer to quote the Extended Source Catalog (XSC) when a match 
is available, but also make use of the Point Source Catalog (PSC). The last two columns in Table 1 
list the IDs assigned to sources by Oyabu et al. (2005) using observations from ISO. Table 2 lists the 
fluxes reported in the SWIRE, 2MASS, and SDSS catalogs along with our own flux measurements 
as described above. 



4.3. Optical spectra 

We collected spectroscopic data on the optical counterparts to the ISO sources from several 
sources, first of which was our library of Keck spectra obtained as part of the original ISO follow-up 
program which had targeted expected counterparts to the ISOPHOT 170/xm sources. Both low- 
and moderate-resolution spectra were taken with the Echellette Spectrograph and Imager (ESI: 
Sheinis et al. 2000) over several observing runs: 2000 Mar 30 - 31 and 2001 Jan 23 - 24 UT for 
the low-resolution data; and 2001 Feb 27 - 28, 2002 Jan 16 - 17, 2002 Feb 16, and 2002 Mar 15 
for the moderate-resolution data. The low-resolution spectra had previously been used for redshift 
identification of putative ISOPHOT 170/xm counterparts and some of these data were published in 



2 Images of the complete sample are available at: http://ifa.hawaii.edu/~bjacobs/LHonlinefigs.pdf 
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Oyabu et al. (2005). The follow-up high-resolution spectra were obtained in order to make accurate 
emission-line flux measurements. These are published here for the first time. In total, we have Keck 
spectra for 19/40 (48%) of our complete sample. An example of our high-resolution Keck spectra 
is shown in Figure 5. In addition to Keck, we supplemented our redshift and flux measurements 
with public SDSS spectra (Adelman-McCarthy et al. 2007) for an additional 9 of our targets. All 
28 of the Keck + SDSS spectra are published on-line in Figures 5.1-5.28 3 . In addition we list 
spectroscopic redshifts for 2 sources for which the spectra themselves are unavailable (see Table 3). 

Ten of our MLPS-160/mi sources do not have optical spectra. Most of these sources were 
either not properly identified in the earlier ISOPHOT-170/Ltm images, or had several radio and 
K-band counterparts within the ISOPHOT-170/t/m beam where the dominant counterpart that had 
previously been targeted with Keck turned out not to be the correct source. For these 10 sources, 
we have been able to determine fairly accurate photometric redshifts, as described below. 

4.4. VLA 1.4 GHz Radio Data 

Deep 1.4 GHz radio continuum images of the LHEX and LHNW fields were obtained using the 
NRAO VLA in the B-configuration in February 2000 and March-April 2001 as part of the AY110 
and AY121 programs. The angular resolution of the data are ~ 5". The achieved sensitivity of the 
LHEX data is la ~ 15 fiJy while the sensitivity for the LHNW field is a factor of 2 worse because 
of a bright (4.2 Jy) continuum source nearby. The photometry is done using the AIPS task SAD, 
which fits a 2D Gaussian to the brightness distribution, and total integrated flux is reported for 
extended sources in Table 2. A more detailed discussion of the radio data is presented by Oyabu 
et al. (2005). 

5. Results 

The complete photometric data set for each of our 40 MIPS-170/iin sources is presented in 
Table 2. We use these data, along with the measured (spectroscopic) and computed (photomet- 
ric) redshifts listed in Table 3, to construct SEDs for each source and to compute total infrared 
luminosities, which are then used to construct the luminosity function for our complete sample. 

5.1. Photometric Redshifts 

The majority of the SWIRE 160^m detections have spectroscopic data, which were used to 
determine their redshifts. For the 10 objects without spectra (noted with "Phot-z" in Table 3), 



3 See: http://ifa.hawaii.edu/~bjacobs/LHonlinefigs.pdf 
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we calculate a photometric redshift using the photometry from SDSS ugriz-b&nds (except for 
J103341. 28+580221.4, see Table 3), 2MASS J H K s -bands and IRAC 3.6 and 4.5/zm. A X 2 template- 
fitting method {Le Phare) was used following the prescription given in Ilbert et al. (2009). This 
method offers an improvement in photo-z accuracy over previous methods, due primarily to im- 
proved calibration using large spectroscopic samples from VLT-VIMOS and Keck-DEIMOS. The 
best fit redshift values and uncertainties are listed in Table 3. 



5.2. Spectral Energy Distributions and Infrared Luminosity 

The photometry and redshifts were used to construct SEDs {vL u ) for each source, which are 
shown in Figure 6 sorted in order of decreasing luminosity. The SEDs are characteristic of what has 
previously been observed for infrared-selected galaxies, with the most luminous sources showing 
a dominant "infrared bump" presumably due to thermal dust emission, and an "optical bump" 
due to thermal emission from stars. Although the mid-infrared sampling is relatively sparse, it is 
also possible to see the effects of emission from polycyclic aromatic hydrocarbons (PAHs) in the 
mid-infrared at A rcst ~ 4 — 12/xm, and silicate absorption at A res t ~ 10//m. 

The SEDs displayed in Figure 6 also show template fits to the data. The MIPS and 8/jm points 
are fit to a library of SED templates by Siebenmorgen & Krugel (2007), and the best fit is shown as 
solid line. The dotted line in Figure 6 represents a stellar evolution model fit to the UV-NIR data 
which is used to estimate stellar masses (see below). To estimate each source's total IR luminosity, 
-^ir(8 — 1000/um), we use the prescription described by Siebenmorgen & Krugel (2007). The use 
of this model of SED fitting to estimate IR luminosity over others, such as Chary & Elbaz (2001) 
or Dale & Helou (2002), is advocated by Symeonidis et al. (2008) largely due to the tendency of 
these models to underestimate the peak of the FIR luminosity as represented by the 160/xm flux. 
This tendency of Siebenmorgen & Krugel (2007) to better fit the FIR data held true for our 160/um 
sample as well. Table 3 lists the computed infrared luminosity for each source. 

5.3. Infrared Galaxy Luminosity Function 

In addition to infrared luminosity, we use each source's flux at 160/t/m, and redshift (Hq = 75 
km s _1 Mpc -1 , J7 m = 0.3 and fl\ = 0.7 to calculate luminosity distance) to calculate a LF 
from our sample. Figure 7 shows the infrared LF resulting from our observations of the LH in 
comparison with the LF in the local Universe, previously determined from the IRAS Revised Bright 
Galaxy Sample (RBGS) all-sky survey, which has median and maximum redshifts: z = 0.008 and 
z = 0.09 (Sanders et al. 2003). The LF density values and uncertainties plotted in Figure 7 are 
listed in Table 4. We divide the data into bins of log(Lm/ Lq) = 0.4 in size, corresponding to 
steps of one in absolute magnitude. Note that SWIRE detected galaxies, J103258. 0+573105 and 
J105349. 60+570708.1 at 160/im, but we do not include them in our analysis because their calculated 
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total infrared luminosities are log(LiR/L ) = 9.55 and 8.81 which results in them falling as the lone 
galaxy in their respective luminosity bins. The volume at which the survey is sensitive to galaxies 
below this luminosity range is small, so we restrict our attention to higher luminosity sources. 

In comparing the LF of our LH data with the RBGS we are comparing a narrow deep survey 
with 40 galaxies to a wide local survey with several hundred galaxies. This means that the relative 
significance of each galaxy is higher for the LH, so it is important to have accurate luminosity 
estimates. On the other hand, the one-magnitude bins have the effect of mitigating uncertainties 
in luminosity, since they can include galaxies with luminosity estimates differing by as much as a 
factor of 2.5. This is particularly significant for the objects in the sample that have luminosity 
estimates using photometric redshifts, since their luminosity is less certain than the rest of the 
sample. The photometric redshifts and their 68% uncertainty limits are noted in Table 3, and these 
limits are used to estimate uncertainties in their luminosities. In addition to photometric redshift 
uncertainties, the errors in matching vL v to a model infrared luminosity become less important 
when the data are binned as described. The space density of galaxies within each bin is calculated 
using the method developed by Schmidt (1968), which accounts for the fact that in a flux-limited 
sample a larger range of luminosities is observable at small distances than at greater distances. He 
proposes using a measure of the volume that a particular flux measurement samples, given the flux 
limits of the survey. For example, a galaxy at redshift z = 0.5 with a flux at 160/xm of 170mJy 
could have been seen at greater redshift (and hence represent a larger volume), since the flux limit 
of the sample is 120mJy. This volume sampling is characterized by the V/V m&x parameter, where 
V is the volume corresponding to the redshift actually observed, and V nmx is the maximum volume 
over which it could be observed. A mean V/V max value of 0.5 within a luminosity bin indicates an 
even distribution of galaxies within the total volume sampled in that bin. 



6. Discussion 

6.1. Evidence For Possible Evolution in the Luminosity Function 

The space density of galaxies with infrared luminosity (8 - 1000/um) in the range log (£ir/ Lq) = 
10 — 12 appears to be consistent between the RBGS and our LH sample. In particular, Sanders et 
al. (2003) fit a broken powerdaw to the RBGS sample. At log (L m /L e ) = 9.5-10.5 the RBGS is fit 
with $(L) oc L-°- 6±ai , and at log (L lR /L Q ) = 10.5-12.5 the powerdaw is: $(L) oc L~ 2 - 2±01 . The 
Lockman Hole data for luminosities log (Lir/L ) < 12, agree within their errors to these power- 
laws. This concurrence is to be expected given the relatively low redshifts sampled in these lower 
luminosity bins. At log (Ljr/Lq) > 12.0 the situation changes. The co-moving space density of 
ultraluminous infrared galaxies (ULIRGs: log (Lir/L ) > 12.0) in the log (Lir/L ) = 12.0 — 12.4 
luminosity bin is ~ 7x higher in the LH than in the RBGS. The median redshift of the ULIRGs in 
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the LH in this luminosity bin is z = 0.51. In the highest luminosity bin, log (Ltr/£©) = 12.4 — 12.8, 
the median redshift of the 4 LH galaxies in this bin is z = 0.71. To compare the co-moving space 
density with the RGBS in this bin we extrapolate the RBGS power-law to log (Ltr/Lq) = 12.6 
and find that the density in the LH sample is ~ 11 x higher. 

Our new results for the LF of the most luminous infrared galaxies in the LH are consistent 
with strong evolution in the co-moving space density of ULIRGs. If we assume pure space-density 
evolution of the form (1 + z) n , our new results for the LH imply n ~ 6 ± 1. This is similar to what 
was found in an earlier study of the infrared luminosity function of ULIRGs by Kim & Sanders 
(1998), where the co- moving space density of ULIRGs in the IRAS 1-Jy sample (mean z ~ 0.15), 
was found to be ~ 2x larger than the local space density of ULIRGs in the RBGS (mean z ~ 0.05), 
implying n = 7.6 ± 3.2 . Our new results are also consistent with a recent determination of the 
extragalactic 250/um luminosity function by Dye et al. (2010), which shows a "smooth increase" 
with redshift of a factor of 3.6 x in the co- moving space density of luminous infrared sources between 
z = and z = 0.2, corresponding to n = 7.1. 

Deeper far-infrared surveys currently underway with Spitzer and Herschel will eventually 
allow us to determine whether the strong evolution observed for the most luminous infrared ex- 
tragalactic sources in the relatively nearby universe continues out to higher redshift. For now, we 
simply note that if we assume similar strong evolution, e.g. (1 + z) 6 , in the ULIRG population out 
to higher redshifts, our results would imply a co-moving space-density of ULIRGs that is ~700x 
larger at z ~ 2 compared to the value at z = 0. Is there evidence for such a large population 
of ULIRGs at high redshift? The answer seems to be yes. There is a population of faint submil- 
limeter sources detected by the Submillimeter Common User Bolometer Array (SCUBA) on the 
James Clerk Maxwell Telescope (JCMT), which has been interpreted variously as exotic objects, or 
ULIRGs at high redshift (Smail et al. 1997; Hughes et al. 1998; Barger et al. 1998; Lilly et al. 1999). 
Lilly et al. argued that these objects are indeed ULIRGs at z ~ 2. Subsequently, Chapman et al. 
(2005) measured a range of spectroscopic redshifts, z = 1.7 — 2.8 for a sample of 73 submillimeter 
galaxies, and suggested an evolution in number density of three orders of magnitude for ULIRGs 
between z = and z ~ 2.5. Our results for ULIRGs in the LH, when extrapolated out to z = 2 — 2.5 
are then consistent with the hypothesis that the SCUBA submillimeter sources are indeed ULIRGs. 



6.2. Galaxy Properties 

To achieve a better understanding of the processes responsible for the observed infrared emis- 
sion and the nature of the galaxies in our MIPS-160/t/m sample, we use our UV-NIR imaging data 
and optical spectra to determine galaxy morphology and masses, and spectral types, respectively. 
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6.2.1. Imaging: Morphology and Masses 

In order to develop a picture of the morphologies, and to gain an indication of the prevalence 
of merging/interacting galaxies in the sample we examine their UV-NIR images. We compile color 
composite images of the sources from those available through the Finding Chart section of the SDSS 
DR7 website, and show these in luminosity order in Figure 8. The zoom on these cutouts is scaled 
so that each box is lOOkpc on a side. At high redshifts and thus high zoom, the image quality of the 
SDSS charts is low, so for sources with z > 0.3 we display stacked g'r'i' images from KPNO when 
available. A brief description of the galaxy morphologies is presented in Table 5. Many of the higher 
luminosity sources with log(LiR./L ) > 11.5 exhibit features suggestive of interactions/mergers, 
such as multiple cores and/or tidal tails. At luminosities lower than log(LiR,/L ) < 11 the large 
majority of sources appear to be mostly unperturbed spirals. These trends are consistent with 
previous studies of local samples of LIRGs and ULIRGs (e.g. Sanders &; Mirabel 1996), which have 
shown that strong interactions and mergers appear responsible for triggering the most luminous 
infrared sources. 

Stellar masses for each of the MIPS-160/mi sources are listed in Table 5. The masses were 
computed by fitting the UV-NIR SEDs using Le Phare (Ilbert et al. 2010) and assuming a Chabrier 
(Chabrier 2003) initial mass function (IMF). The mass range is log (M/M ) ~ 10.0 - 11.5 cor- 
responding to ~ 0.5 — 3M*. Higher mass systems are more likely to be associated with higher 
infrared luminosity. 

6.2.2. Spectroscopy: Extinction, Abundances and Spectral Types 

Our high- resolution Keck/ESI spectra, supplemented by SDSS spectra and four low-resolution 
ESI spectra, allow us to measure robust spectral types for 25 of the 160/im sources in our sample. 
An example of these spectra was shown in Figure 5. Spectra for the 25 sources with spectral types 
as well as those from three sources for which we have data but were unable to measure spectral 
types are available in the online edition of the Journal. After carefully accounting for the effects 
of stellar absorption and applying an extinction correction (median E(B — V) = 0.7), we classify 
the spectra as H II-region-like, or star- forming (H); composite, or star formation + an AGN (C); 
Seyfert (S); or LINER (L), based on the classification scheme proposed by Kewley et al. (2006) (see 
Figure 9). Table 3 lists the measured spectral types and extinctions for each galaxy when they are 
available. 

Dividing the subsample with spectral types into luminosity bins, we find that 15 of 19, (79%), 
of galaxies in the log(LiR,/L ) < 11 bin have H II-region-like spectral type, consistent with star 
formation as the dominant source of excitation. This is close to the H II-region-like fraction of 
nearby galaxies selected at 60/mi (~70% - Veilleux et al. 1995; Yuan et al. 2010). The other four 
galaxies include one Seyfert and one LINER and 2 objects with mixed types that suggest a "com- 
posite" starburst-AGN mixture of excitation. Only six galaxies with high resolution spectroscopy 
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(and hence derived spectral types) have log(LiR,/L ) > 11. Three (3/6 = 50%) have H II-region- 
like spectral type, one is a Seyfert, one is a Seyfert/LINER, and one is a "composite" mixture of 
starburst /LINER excitation. Although the fraction of galaxies with H II-region-like spectral types 
decreases at higher infrared luminosity (similar to what is observed for nearby galaxies selected at 
60/xm), the number statistics in this high luminosity bin are too low to draw conclusions about 
the fractions of different spectral types. In an attempt to provide additional information on the 
spectral types of our high infrared luminosity sources, we have employed a new technique developed 
from studies of SDSS galaxies (Smolcic et al. 2008) that maps UV/optical continuum colors onto 
the spectral line diagnostic diagram. This method is described in the Appendix, where the SDSS 
photometry for all of our MIPS-160^m sources is used to derive "P1,P2" photometric spectral 
types for each source, following the prescription given by Smolcic et al. (2008). These results both 
confirm the large H II-region-like fraction among the lower luminosity infrared sources, and show 
that composite and AGN spectral types appear to increase among the highest luminosity sources. 

Because our spectra also contain the [O II] AA3727, 3729 doublet, we are able to estimate 
gas-phase oxygen abundances for these systems. Where available, these are listed in Table 5, using 
the robust [N II] /[O II] diagnostic of Kewley et al. (2002). To put these in context, we also used the 
measured K s data to compare to the luminosity- metallicity relation in the NIR (Salzer et al. 2005). 
For the 8 systems that have sufficient information (upper-branch R23 gas abundances and measured 
luminosities; see Rupke et al. (2008) for more on the methodology), we find that 5 follow the L — Z 
relation of normal galaxies. Three others have higher luminosities than the data threshold, and 
appear to be slightly below the L — Z relation (by 0.1 — 0.2 dex), as found for other infrared-selected 
objects at high luminosity (Rupke et al. 2008). 

6.2.3. Radio-FIR Correlation 

The measured 1.4 GHz radio continuum fluxes are converted to 1.4 GHz radio power Li.4GHz 
assuming a spectral index of a = +0.75 4 , and they are plotted as a function of redshift on the 
left panel of Figure 10. The observed 1.4 GHz radio power range between 10 20 and 10 24 3 W 
Hz -1 (see Table 3), similar to the IR-selected galaxies in the local universe studied by Yun et al. 
(2001b), and none of the sources has sufficient radio power to be classified as a "radio-loud" object. 
The most luminous infrared sources also tend to be those with the largest radio luminosities, i.e. 
log(Li.4QHz) = 23.0 — 24.4, equivalent to the radio powers typically seen among Seyfert galaxies, 
and thus the presence of a low luminosity AGN cannot be ruled out based on these radio powers 
alone. 

The well-known correlation between the measured infrared luminosity and radio power for star 



4 Spectral index a is defined as S v = Sq{^-) 
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forming galaxies is often quantified using the ratio commonly referred to as "g- value" 

q = log [(FIR/3.75 x 10 12 H Z )/5 1 . 4G Hz] (1) 

where FIR is the far-infrared flux density and Si.4GHz is m W m~ 2 Hz -1 (Condon 1992; Yun et al. 
2001b). We computed these q- values using the Li.4GHz derived above and I/fir computed from the 
best-fit SED models integrated between A = 40 and 500 |im, where the wavelength range has been 
chosen to match the original definition of Lfir used to compute u q" . As shown on the right panel 
in Figure 10, the derived g-values of the LH 160 fim sources fall between 1.6 and 3.0, suggesting 
that most of these sources follow the same radio-FIR correlation as the local star forming galaxy 
populations. Some of the low redshift (z < 0.15) sources appear to have q- values on the high end 
of the local population. These are also the sources with the largest angular size, and the VLA 
measurements are likely under-estimates as a consequence. None of the LH 160 /im sources has a 
q- value less than 1.6 and thus a clear evidence for a radio-loud AGN. 

7. Summary 

We have made use of Spitzer-SWIRE imaging data of two ~ 0.5deg 2 fields (LHEX and LHNW) 
in the Lockman Hole, to identify a complete sample of 40 MIPS-160^m selected extragalactic 
sources, with S*i6o > 120mJy. In combination with Keck spectroscopy and photometry from SDSS 
and 2MASS, we have obtained redshifts and infrared luminosities, and have attempted to charac- 
terize the host galaxy properties for all of the objects in the sample. The luminosity function for 
the MIPS-160 / um sample has been compared with the "local" (z < 0.05) luminosity function of 
FIR galaxies previously derived using the IRAS all-sky survey. 

Our main results can be summarized as follows: 

(1) The complete 5i6o > 120mJy sample contains 40 galaxies with infrared luminosities in the 
range log(LiR,/L ) = 8.81 — 12.74, with a maximum redshift: z = 0.80. The cumulative source 
counts down to 120mJy are estimated to be 1.2 x 10 5 sources sr _1 at 160^tm. 

(2) The luminosity function of the sources with log(Lni/ 'Lq) ~ 9.5 — 11.5 is similar to that 
found previously for infrared galaxies in the IRAS 60/im local galaxy sample. 

(3) The co- moving space density of the MIPS-160/zm galaxies with log(Lm/ L Q ) > 12 is ~ 10 x 
higher than that for the local infrared galaxies with similar infrared luminosities found in the 
IRAS RBGS. Assuming pure number density evolution proportional to (1 + z) n , these results give 
n = 6± 1, which implies strong evolution of the most luminous infrared sources, in contrast to little 
or no evolution observed in the number density of lower luminosity objects with log(LiR/L ) < 11. 

(4) The host galaxy masses for our sample are in the range log(M*/MQ) ~ 10.0 — 11.5 
(0.5 — 3M*), with evidence for an increase in host mass from a mean of 10.7 to 11.0 for ob- 
jects with infrared luminosities below and above log(LiR,/L ) = 11, respectively. 

(5) The morphology and spectral types for our flux-limited sample of 160/um-selected sources 
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generally agree with what has been observed locally for 60^m-selected samples. At log(LiR,/L ) > 
11, the fraction of disturbed and/or merger systems and the fraction of objects with "composite" 
and/or AGN spectral types increases with increasing Ljr. At log(LrR/L0) < 11 most objects ap- 
pear to be either unperturbed spirals and/or weakly interacting systems with spectral types typical 
of HII regions. 

(6) None of the LH sources has sufficient radio power to be classified as a "radio-loud" object. 
However, the most luminous infrared sources also tend to have the highest radio luminosities, i.e. 
log(Li.4GHz) = 23.0 — 24.4, equivalent to the radio powers typically seen among Seyfert galaxies, 
and thus the presence of a low luminosity AGN cannot be ruled out based on these radio powers 
alone. 
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A. Using rest-frame colors to assign galaxy types 

In an attempt to determine the spectral type of the 12 galaxies in our sample for which we lack 
optical spectra, we make use of a photometric (rest-frame color based) method extensively studied in 
Smolcic et al. (2008). Smolcic et al. (2006) have shown that principal component rest-frame colors 
(PI and P2 hereafter) drawn from the modified Stromgren photometric system (3500-5800 A; 
Odell et al. 2002) essentially trace the position of a galaxy in the BPT diagram (Baldwin, Phillips 
k, Terlevich 1981; Kewley et al. 2001, 2006; Kauffmann et al. 2003). Hence, they can be used as 
efficient tracers of galaxy type, such as low-luminosity AGN (Seyfert and LINER), star forming, 
and composite galaxies (see Smolcic et al. 2006, 2008 for more details). Here we extend the method 
developed by Smolcic et al. (2008), that utilizes only the PI color to disentangle star forming from 
AGN galaxies, to a probabilistic approach that uses both, PI and P2 colors, and we adapt it to an 
IR- (rather than radio-) selected sample. 

A.l. Derivation of PI, P2 colors 

Given that we have SDSS ugriz photometry for the complete MIPS-160/im sample, we derive 
the PI and P2 colors for all 40 galaxies in our sample by fitting their SEDs (encompassed by the 
SDSS ugriz photometry) with 100,000 spectra from the Bruzual & Chariot (2003) stellar population 
synthesis model library. Before performing the x 2 minimization fit we redshift all the model spectra 
to the galaxy's spectroscopic redshift. The colors are then computed from the best fit model 
spectrum (see Sec. 4.2. in Smolcic et al. 2008 for more details about the SED fitting). 

To assess the accuracy of photometrically synthesized colors for IR-selected galaxies, we have 
derived the (P1,P2) colors via SED fitting (in the same way as described above) for an IR-selected 
control sample with available (P1,P2) colors computed independently from their spectra. The 
control sample, limited to a redshift range of 0.04 to 0.3, contains ~ 1350 galaxies drawn from 
the SDSS DR1 "main" spectroscopic galaxy sample matched to the IRAS Faint Source Catalog 
(see Obric et al. 2006 for details about the cross-correlation of the catalogs). The spectroscopically 
derived (P1,P2) colors have been synthesized by Smolcic et al. (2006) by convolving the SDSS high- 
resolution spectra with the Stromgren filter system. They estimated that these spectroscopically 
derived rest-frame colors are accurate to 0.03 mag. 

In Figure 11 we show the difference between the photometrically and spectroscopically derived 
(P1,P2) colors for our IR-selected control galaxies as a function of the photometrically derived PI 
color. We use the median offset as a function of PI to correct the photometrically derived colors, i.e. 
to scale these to the SDSS spectroscopic system. The distribution of the (P1,P2) color differences 
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after the corrections have been applied are shown in Figure 12. As expected, the corrections have 
eliminated systematic effects. Furthermore, as the accuracy of the spectroscopically derived colors 
has been shown to be 0.03 mag (Smolcic et al. 2006), the accuracy of the photometrically derived 
(P1,P2) colors is likely better than 0.12 and 0.03, respectively. It is remarkable that the accuracy 
of the P2 color derived via SED fitting is comparable to that of the spectroscopically derived color. 

A. 2. Classifying the Lockman galaxies using the P1-P2 color method 

In Figure 13 we show the PI vs. P2 the distribution for our 40 galaxies in our MIPS-160/xm 
sample (their P1,P2 colors were corrected for systematic offsets as described in the previous section; 
filled dots). To assess the nature of these 40 galaxies for each one we compute the probability that 
it is a star forming, composite, AGN or absorption galaxy given its (P1,P2) rest-frame colors. 
The probability is computed based on the underlying distribution of our ~ 1300 control (SDSS- 
IRAS; 0.04 < z < 0.3) galaxies (spectroscopically divided into absorption, star forming, AGN and 
composite galaxies using the standard diagnostics; Baldwin, Phillips & Terlevich 1981; Kewley 
et al. 2001; Kauffmann et al. 2003) in the P1-P2 plane as follows. We bin the P1-P2 plane in 2 
dimensions. The size of the bins is taken to be about 2 times the photometric color uncertainty (see 
previous section). For each (P1,P2) bin we then calculate the probability as the ratio of the number 
of each (spectroscopically classified) galaxy type relative to the total number of control galaxies 
in that bin. Given the photometrically synthesized (P1,P2) colors for our Lockman galaxies we 
can then access the probability of each galaxy being a star forming, AGN, composite or absorption 
galaxy. These probability contours are shown in Figure 13, and our results are summarized in Table 
6. 
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The lochirian Hole Field 
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Fig. 1. — SWIRE and ISO coverage of the Lockman Hole. The black box denotes the area covered 
by SWIRE MIPS observations, the magenta box denotes the SWIRE IRAC channels 1 and 3, 
and the cyan box denotes SWIRE IRAC channels 2 and 4. The red square marks the location of 
observations of the ISO LHEX region, the red polygon towards the right of the image marks the 
ISO LHNW region. 
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Fig. 2. — Left panel: Distribution of fluxes of the ~500 sources in the SWIRE LH 160/um catalog. 
Each bin is lOmJy wide. Right panel: Integral number counts versus flux of sources in the SWIRE 
LH 160^m catalog (horizontal axis ranges from 50 to 300 mJy). Vertical bars denote Poisson 
uncertainty. 
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Fig. 3. — Left panel: LHEX subregion of the LH. The coverage of LHEX is roughly: a = 
10 h 49 m 30 s - 10 h 54 m 30 s and 5 = 57°00'00" - 57°44'00". The circles denote the location of the 
160/im sources. Right panel: LHNW subregion of the LH denoted by the large polygon. The 
coverage of LHNW is roughly: a = 10 h 31 m 00 s - 10 h 37 m 00 s and 8 = 57°20'00" - 58°10'00". 
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Fig. 4.— Images of source J103320. 32+574913.6 at g', r' , i'-bands, and 3.6, 4.5, 5.8, 8.0, 24, 
70, 160/um. The red, blue, and green circles indicate the position of the centroid of each de- 
tection in the 160/um, 70/rni, and bandmerged IRAC+24 / um catalogs, respectively. The radii 
of the circles are set to 40", 20", and 5" in all images (for all sources), and are meant to 
be rough indicators of the size of the PSF for the 160/mi, 70//m, and 24^m sources. Figures 
4.1 - 4.40 display cutout images for the entire sample of 160/im sources and are available at: 
htt p : / /if a . hawaii . edu/ ~ b j acobs /LHonlinefigs . p df 
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Fig. 5. — Example spectrum of J103320. 32+574913. 6 (see Table 3 for the source of the data). The 
wavelength range is 3650 - 6800A, and the vertical lines mark the positions of [OII]3726/3729, Call 
H & K, H/3, [OIII]5007, Nal D, Ha, and [NII]6583. Figures 5.1- 5.28 display spectra for the 28 
sources for which we have data, see: http://ifa.hawaii.edu/~bjacobs/LHonlinefigs.pdf 



-22- 



ion 
10'* 

1C S 
10 ,a 
10" 1 

l0 ld 

■■■J 

10* 
101* 

f> 10' 1 
/ 10'* 
10 s 

10" 
10 11 
lO 1 * 

10* 

in* 

10 13 

10" 1 
10 io 

10 s 
10* 



(12-ZO) 




!■ ■ ■ ■ — \ — 



(11.40 




I r 




Jinjr_ijB.8*5TSS4* 
(fig) 



... 



( I I -37} 




:-:s.H'); 



J M34I ,ZBiiHG221 4 



(IZ.Z&) 




— 1, , , , — |„ — - 



(n,9&) 



■1105056.80+571*31,2 




j:>M4;' .40+5724*4,7 

P1.H) 



4>v « 



(10,87) 




jicjjii.ttiti/jna 
f\ {12.4*) 



i >■ 



{12.12^ 



Y — I — 



[11.04) 



/ 




J10M4 .9atS707ft6.7 
(10.7$) 



V 

ft 



100 10 1 ioo 10 i 100 io 

Pest Wavelength (^rn) 



100 10 1 



Fig. 6. — SEDs of SWIRE 160/im sources in vL v in units of solar luminosity versus rest wavelength (increasing to the left). 
The solid line represents a Sicbcnmorgcn & Krugcl (2007) model fit to the MIPS and 8/im data, shown as diamonds. The dotted 
line represents the stellar component of the luminosity and is fit to the shorter wavelength data (Ilbcrt et al. 2010). Sources 
are sorted by decreasing IR luminosity. We list the total-IR luminosity in each panel in units of log(.LiR,/Z/0) in parentheses 
below the source ID. 
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Figure 6 (Continued). 
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Fig. 7. — Infrared galactic LF from observations of the LH by SWIRE are shown as filled diamonds 
with 68% uncertainty bars. The RBGS data are shown as pluses fit by two powerdaws (Sanders et 
al. 2003). At log(L m /L ) = 9.5- 10.5 the RBGS is fit with oc L -0 - 6 * - 1 , and at log(L m /L ) = 
10.5 - 12.5 the powerdaw is: $(L) oc L _2 - 2±01 . 
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Fig. 8. Stacked images of the LH sources sorted by luminosity. Most images are from SDSS finding charts, however for 

sources at z > 0.3 we show stacked KPNO g'+r'+i' images when available, because these reveal more detail than SDSS (these 
images are marked with an asterisk). Images are oriented with North up and East to the left, with dimensions of 100 kpc X 
100 kpc. We list the total-IR luminosity for each source in units of log(LiR/LQ) in parentheses following the source ID. 
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Figure 8 (Continued). 
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Fig. 9. — Optical spectral diagnostics, following the formalism of Kewley et al. (2006), for those 
sources with optical spectra and for which we were able to measure the relevant diagnostic lines. 
The dotted lines divide the diagrams into spectral classification types, and the solid line denotes 
the limit to pure star formation. 
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Fig. 10. — Left panel: 1.4 GHz radio power of the 34 sources with available radio data shown 
as a function of their redshifts. Right panel: The radio-FIR correlation (/-value of the same 34 
sources are shown. The solid horizontal line corresponds to q = 2.34, which is the mean of the local 
star forming galaxies, while the lower and upper dashed lines mark the radio-excess and IR-excess 
objects, respectively (see: Yun et al. 2001b). 
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Fig. 11. — Offsets between the photometrically and spectroscopically derived PI and P2 colors as 
a function of the former for a sample of ~ 1350 SDSS-IRAS galaxies (small gray dots). The former 
were obtained via SED fitting to the SDSS ugriz photometry using 100,000 model spectra from 
the Bruzual & Chariot (2003) library. The latter were derived from SDSS spectra as described in 
Smolcic et al. (2006). The large dots represent median offsets as a function of the photometrically 
derived PI color. The error bars show the interquartile ranges. 
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Fig. 12. — Distribution of the photometrically and spectroscopically derived PI and P2 colors after 
the systematic offset correction (shown in Figure 11) was applied. The number of objects, mean, 
and standard deviation, as well as the best Gaussian fit are also shown in the panels. 
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Fig. 13. — PI vs. P2 color plane. The colors for the Lockman galaxies (after the corrections for 
systematics have been applied) are shown by large dots. Sources with 11 < logLiR/L < 12 
and logLift/L© > 12 are encircled with dashed and full lines, respectively. Contours represent 
the absorption (i.e. galaxies with no emission lines; NOEML), star forming, composite and AGN 
probability levels given a (P1,P2) combination (see text for details). The uncertainties in the PI 
and P2 colors obtained via SED fitting are shown in the top left of the panel (see also Figure 12. 
Note, however, that the P2 error is likely lower than indicated in the panel (see text for details). 
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Table 3. Redshift and Infrared Luminosity 



# 




Source 




z 


z-Ref. 


log(Li R ) 
(Lq) 


log(Ll.4GHz) 

(W Hz" 1 ) 


1 


SWIRE3. 


J103237.44+580845.9* 





.42 + 0.02 


Phot-z 


12.28 + 0.05 




2 


SWIRE2. 


.24_J 103249.4+573707 





.115 


ESI/ech 


10.67 


22.40 


3 


SWIRE3 


J103253.94+580633.0 





.073 


SDSS 


10.17 




4 


SWIRE2. 


.24.J103258.0+573105 





.046 


SDSS 


9.55 


< 20.60 


5 


SWIRE2. 


_24_J103314.8+573110* 





_i_n 99 

•8o±8:?l 


Phot-z 


_i_n 9K 
12.48l°' 2 « 


< 23.77 


6 


SWIRE3. 


.J103320.32+574913.6 





.240 


ESI/ech 


11.37 


23.02 


7 


SWIRE3. 


.J103327.90+574534.4 





.152 


SDSS 


10.60 


< 21.78 


8 


SWIRE3 


J103341. 28+580221.4* 





-i_n or 
oz -0.08 


Phot-z 2 


1 2 62 + 


24.05 


9 


SWIRE3. 


.J103358.73+574317.1 





.263 


ESI/ech 


11.46 


23.49 


10 


SWIRE2 


24J103358.9+572952 





.469 


ESI/ech 


12.04 


24.24 


11 


SWIRE3. 


.J103515.83+573337.4* 





.16 + 0.02 


Phot-z 


10.87 + 0.12 


< 21.84 


12 


SWIRE3. 


.J103526.79+575147.6 





.595 


ESI/low-d 


12.25 


23.85 


13 


SWIRE2. 


_24_J103538.8+573546* 





.69 + 0.03 


Phot-z 


12.66lg;°5 


24.36 


14 


SWIRE3. 


J103539.24+574243.9 





.113 


ESI/ech 


10.58 


21.93 


15 


SWIRE2 


24 J103557. 1+572234* 





4K+0-02 
4D -0.26 


Phot-z 


1 1 qo+0.04 
11.9Z_ g4 




16 


SWIRE3. 


J103603.97+574812.5 1 





.511 


B 1 


12.20 




17 


SWIRE3 


J103606.48+574702.4 





.044 


SDSS 


9.95 




18 


SWIRE3. 


.J103653.50+575442.5 





.102 


SDSS 


10.46 




19 


SWIRE3. 


.J104948.86+573458.2 





.135 


SDSS+ESI/low-d 


10.89 


22.20 


20 


SWIRE3. 


.J104956.06+571440.4* 





.07 + 0.03 


Phot-z 


q qa+0.37 


< 20.87 


21 


SWIRE3 


J105041.96+570706.7 





.091 


ESI/ech 


10.76 


22.18 


22 


SWIRE3. 


.J105052.41+573506.9 





.028 


SDSS+ESI/ech 


9.64 


20.64 


23 


SWIRE3. 


.J105056.60+571631.2 





.397 


ESI/low-d 


11.69 


23.55 


24 


SWIRE3 


J105100.42+574114.9 





.071 


ESI/ech 


10.12 


< 21.12 


25 


SWIRE3. 


J105113.41+571425.9 





.362 


ESI/ech 


11.90 


23.76 


26 


SWIRE3. 


J105143.75+572936.9 





.081 


ESI/ech 


10.18 


22.01 


97 


o VV J-TVCjO. 




u 


1 1R 
■ lOU 


OUDj 


1U.OU 


91 ^7 


28 


SWIRE3. 


J105151.64+570935.7 





.608 


ESI/low-d 


12.12 


23.44 


29 


SWIRE3 


J105207. 16+570745.5 





.123 


ESI/ech 


10.73 


22.05 


30 


SWIRE3. 


.J105225.75+570153.6 





.061 


ESI/ech 


10.16 


21.66 


31 


SWIRE3. 


.J105242.40+572444.7* 





lq +o.oi 


Phot-z 




22.57 


32 


SWIRE3. 


.J105252.76+570753.7 





.163 


SDSS+ESI/low-d 


10.99 


22.53 


33 


SWIRE3 


J105256.85+570825.6 





.080 


ESI/low-d 


10.40 


21.90 


34 


SWIRE3. 


J105301. 36+570543.1 





.080 


ESI/ech 


10.64 


22.14 


35 


SWIRE3. 


.J105308.32+570645.6* 





7 g+0.14 
,o -0.26 


Phot-z 


12 74 +o.i9 


< 23.86 


36 


SWIRE3. 


.J105314.84+574137.6* 





91 +0.05 
z± -0.03 


Phot-z 


11 04+ ' 21 


22.41 


37 


SWIRE3 


J105318.94+572140.7 





.133 


SDSS 


10.70 


22.34 


38 


SWIRE3. 


.J105320.92+571433.2 





.080 


SDSS 


10.03 


< 20.89 


39 


SWIRE3. 


.J105349.60+570708.1 1 





.006 


NED 1 


8.81 


19.84 


40 


SWIRE3. 


J105432.31+570932.4 





.068 


SDSS 


10.44 


21.86 



Note. — Sources in the LHNW and LHEX fields listed in order of increasing right ascension, with 
columns giving a counter number, source ID, redshift, infrared luminosity (8-1000/im) in units of Lq, and 
radio power Li 4qhz m units: log(W Hz -1 ). The sources of the spectra listed in the z-Ref. column 
are: from ESI on Keck II in low-dispersion (low-d) and echellette (ech) mode, from Amy Barger's private 
communication (2003) (B), from SDSS, and from NED. 
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* Sources without existing optical spectra, for which we have determined photometric rcdshifts (see text 
and Figure 6). The computed phot-z value, and 68% uncertainty is listed. This uncertainty was propagated 
when computing the listed infrared luminosity. 

1 Thirty of the forty sources in our sample have spectroscopic redshifts, and 28 of these have spectra, which 
are displayed in Figures 5.1 - 5.28. The spectra for sources J103603.97+574812.5, and J105349. 60+570708.1 
are not available. 

2 For the case of J103341. 28+580221. 4, we prefer the SWIRE optical data due to confusion with a fore- 
ground spiral and use these rather than the SDSS data to calculate its photometric redshift. 
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Table 4. Luminosity Function 



log(L IR ) 




<7<J> 


v/v max 


°V/V max 


Number 


Median z 


(A=>) 


(Mpc -3 mag -1 ) 












9.8 


3.1 x 10 -3 


1.8 x 10 -3 


0.53 


0.11 


3 


0.044 


10.2 


1.4 x 10 -3 


5.6 x 10~ 4 


0.64 


0.09 


6 


0.077 


10.6 


7.2 x 10~ 4 


2.3 x 10~ 4 


0.59 


0.07 


10 


0.114 


11.0 


8.3 x 10~ 5 


3.7 x 10~ 5 


0.68 


0.11 


5 


0.163 


11.4 


8.7 x 10~ 6 


6.1 x 10~ 6 


0.66 


0.17 


2 


0.251 


11.8 


2.4 x 10~ 6 


1.4 x 10~ 6 


0.65 


0.13 


3 


0.397 


12.2 


1.4 x 10~ 6 


6.2 x 10" 7 


0.57 


0.09 


5 


0.511 


12.6 


3.6 x 10~ 7 


1.8 x 10~ 7 


0.58 


0.10 


4 


0.725 



Note. - The galaxy sample is placed into luminosity bins of 0.4 in the range 
log(LiR/L ). We list the space density <!>, the uncertainty in this value, volume sam- 
pling parameter V/V mSLX , uncertainty in V/V max , the number of galaxies in each bin, 
and the median redshift. 
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Table 6. Spectral Type Probabilities from P1P2 Analysis 



# Source log(Lin) PI Color P2 Color SF Comp. AGN No EML 



(Lq) Probabilities 



35 


SWIRE3. 


.J105308.32+570645.6 


12 74+ - 19 


13 


SWIRE2. 


_24_J103538.8+573546* 


12 66+ ' 05 


8 


SWIRE3. 


J103341. 28+580221.4* 


12 62+ - 10 
-0.15 


5 


SWIRE2. 


.24.J103314.8+573110 


12.48l°-j6 


1 


SWIRE3. 


J103237.44+580845.9 


12.28 + 0.05 


12 


SWIRE3. 


J103526.79+575147.6 


12.25 


16 


SWIRE3. 


J103603.97+574812.5 


12.20 


28 


SWIRE3. 


.J105151.64+570935.7* 


12.12 


10 


SWIRE2. 


_24_J103358.9+572952 


12.04 


15 


SWIRE2. 


24J103557. 1+572234 


11 o 2 + 004 


25 


SWIRE3. 


.J105113.41+571425.9 


11.90 


23 


SWIRE3. 


J105056.60+571631.2 


11.69 


9 


SWIRE3. 


J103358.73+574317.1 


11.46 


6 


SWIRE3. 


.J103320.32+574913.6 


11.37 


31 


SWIRE3. 


J105242.40+572444.7 


11.05t°° B 


36 


SWIRE3. 


.J105314.84+574137.6 


11.04±g;?J 


32 


SWIRE3. 


J105252. 76+570753. 7 


10.99 


19 


SWIRE3. 


J104948.86+573458.2 


10.89 


11 


SWIRE3. 


J103515. 83+573337.4 


10.87+0.12 


21 


SWIRE3. 


.J105041.96+570706.7 


10.76 


29 


SWIRE3. 


.J105207.16+570745.5 


10.73 


37 


SWIRE3. 


J105318.94+572140.7 


10.70 


2 


SWIRE2. 


.24.J103249.4+573707 


10.67 


34 


SWIRE3. 


.J105301.36+570543.1 


10.64 


7 


SWIRE3. 


J103327.90+574534.4 


10.60 


14 


SWIRE3. 


.J103539.24+574243.9 


10.58 


27 


SWIRE3. 


J105150.50+573905.7 


10.56 


18 


SWIRE3. 


.J103653.50+575442.5 


10.46 


40 


SWIRE3. 


.J105432.31+570932.4 


10.44 


33 


SWIRE3. 


J105256.85+570825.6 


10.40 


26 


SWIRE3. 


J105143.75+572936.9 


10.18 


3 


SWIRE3. 


J103253.94+580633.0 


10.17 


30 


SWIRE3. 


J105225. 75+570153. 6 


10.16 


24 


SWIRE3. 


.J105100.42+574114.9 


10.12 


38 


SWIRE3. 


J105320.92+571433.2 


10.03 


20 


SWIRE3. 


J104956.06+571440.4 


o og+0-37 


17 


SWIRE3. 


.J103606.48+574702.4 


9.95 


22 


SWIRE3. 


J105052.41+573506.9 


9.64 


4 


SWIRE2. 


_24J103258.0+573105 


9.55 


39 


SWIRE3. 


.J105349.60+570708.1 


8.81 



0.37 


0.06 


0.37 


0.51 


0.06 


0.06 




* 




* 


* 




* 


* 


* 


* 


* 




0.35 


0.00 


0.08 


0.19 


0.29 


0.44 


0.11 


0.03 


0.54 


0.31 


0.10 


0.04 


0.05 


0.06 


0.43 


0.53 


0.04 


0.00 


-0.51 


0.01 


0.90 


0.03 


0.02 


0.05 




* 




* 






0.33 


0.09 


0.37 


0.51 


0.06 


0.06 


0.22 


0.05 


0.37 


0.51 


0.06 


0.06 


0.49 


0.01 


0.06 


0.33 


0.36 


0.24 


-0.10 


0.01 


0.72 


0.25 


0.02 


0.00 


0.11 


0.07 


0.43 


0.53 


0.04 


0.00 


-0.08 


0.03 


0.72 


0.25 


0.02 


0.00 


-0.38 


0.02 


0.90 


0.03 


0.02 


0.05 


-0.19 


0.02 


0.72 


0.25 


0.02 


0.00 


0.11 


0.05 


0.54 


0.31 


0.10 


0.04 


0.37 


0.07 


0.37 


0.51 


0.06 


0.06 


0.10 


-0.02 


0.34 


0.29 


0.24 


0.12 


0.35 


0.06 


0.37 


0.51 


0.06 


0.06 


-0.12 


0.01 


0.72 


0.25 


0.02 


0.00 


-u. ±u 


u.uo 






n 09 


n no 


0.04 


0.03 


0.54 


0.31 


0.10 


0.04 


-0.05 


0.03 


0.54 


0.31 


0.10 


0.04 


-0.05 


0.02 


0.54 


0.31 


0.10 


0.04 


0.16 


0.04 


0.54 


0.31 


0.10 


0.04 


0.15 


0.04 


0.54 


0.31 


0.10 


0.04 


-0.19 


0.03 


0.72 


0.25 


0.02 


0.00 


0.20 


0.03 


0.54 


0.31 


0.10 


0.04 


0.32 


0.05 


0.22 


0.47 


0.19 


0.12 


0.27 


0.05 


0.37 


0.51 


0.06 


0.06 


-0.27 


0.02 


0.72 


0.25 


0.02 


0.00 


-0.44 


0.01 


0.90 


0.03 


0.02 


0.05 


-0.14 


0.01 


0.72 


0.25 


0.02 


0.00 


0.45 


0.02 


0.22 


0.47 


0.19 


0.12 


0.16 


0.03 


0.54 


0.31 


0.10 


0.04 


-0.30 


0.03 


0.72 


0.25 


0.02 


0.00 


0.07 


0.04 


0.54 


0.31 


0.10 


0.04 


0.29 


0.03 


0.22 


0.47 


0.19 


0.12 


0.42 


0.03 


0.22 


0.47 


0.19 


0.12 



Note. — Sources in the LHNW and LHEX fields listed in order of decreasing Lir. Columns give RA ordered 
counter number, source ID, infrared luminosity (8-1000/im) in units of log(L0), corrected PI color, corrected P2 color, 
Star-Forming probability, Composite probability, AGN probability, and No Emission-Line probability calculated by 
comparing PI and P2 colors with those of the IRAS sample (see Appendix). 
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*Sources J103538. 8+573546 and J103341. 28+580221. 4 suffer from contamination due to partially overlapping fore- 
ground galaxies, which prohibits determining accurate P1,P2 colors. Source J105151. 64+570935. 7 is spcctroscopically 
classified as a Seyfert and is suspected to harbor a QSO, thus the use of P1,P2 color to determine spectral type is 
inappropriate. 



